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OF WING-BODY-TAIL COMBINATIONS 
By George H. Holdaway and William A. Mersman 

SUMMARY 


The techniques of the computing procedure of NACA RM A 53 HI 7 have been 
significantly improved by a new procedure of harmonic analysis using 
Tchebichef polynomials . This Improved method is described in detail with 
illustrations of its two main advantages; these are, simplification of 
the computing procedures, and the provision for a comprehensive check, 
solution which includes a direct check of how well the number of harmonics 
used represent the area-distribution curve. Far the present, no specific 
recommendation can be .made as to the number of harmonics which should be 
used for all configurations in making wave-drag computations; however, 
certain guides Eire given in the concluding remarks of the report. The 
new procedure is also evaluated by camparisons with analytical solutions, 
results from the previous method, and experimental resxzlts. 


IHTRODUCTION 


The computing method of reference 1 has been effectively used to 
estimate the zero-lift drag -rise coefficients of various relatively smooth 
wing -body- tail ccxnblnations (refs. 2, 5, and i)-) . However this original 
method Involves several operations, and the checking procedures, such as 
were used in reference 5^ ars time-consuming and check back to only the 
slopes of the area-distribution ciirves and not to the area-distribution 
curves themselves. 

It is the purpose of this paper to present and to analyze another 
method for representing the slope of area-distribution curves, which will 
allow for a more rapid computation of wave drag and will permit the use 
of an improved method of checking the computations . The basic method for 
computing the wave drag is fundamentally the same as in reference 1 and 
is based on the theory of reference 5* main difference is that 
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Tcheblchef polynomials (refs. 6 and also spelled Chebysiiev) are used 
to represent the Fourier sine series defining the slopes of the area- 
distribution curves . 

The new method is evaluated by comparison of results with analytical 
solutions, results from the previous method, and experimental results for 
Mach numbers up to 1.8. The configurations S'elected, for which experi- 
mental data were available, included models of a trlangilLar-wing 
interceptor-type airplane, a swept-wing interceptor-type airplane, and a 
body- tall configuration with a scoop-inlet duct. 


SiaiBOLS 


An 

CD. 




c 

M.A.C. 

< 3-0 
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M 
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<1 
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coefficients defining the magnitude of the harmonics of a 
Fourier sine series 


zero-lift drag coefficient. 


drag at zero lift 


zero-lift wave-drag coefficient. 


zero-lift drag-rise coefficient. 


qSv 

theoretical wave drag 
at zero lift 

zero-lift drag rise above 
subsonic level 


local chord measured parallel to plane of symmetry 

mean aerodynamic chord of the total wing 

maximum body diameter 

fuselage or body length 

free -stream Mach number 

duct mass-flow ratio ‘ 

mmiber of terms or harmonics used in the Fourier sine series 
a harmonic of the Fourier sine series 
free-stream dynamic press-ure 

projection of Sg on a plane perpendicular to x axis 
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Ss 

S*(x) 

X 


x,y,z 


0 


I 

9 


area formed by cutting configurations with, planes perpendicular 
or oblique to the x axis 

slope of S curves as a function of x, 

' dx 

total wing area 

distance along the x axis measured from the mid-length 
position 

Cartesian coord in ates as conventional body axes 

angle between the z axis and the intersection of the cutting 
planes X with the yz plane 

(See ref. 1 for descriptive sketches and detailed definitions.) 

distance along the x axis measured ftom the mid-length position 
divided by one half of the body length 

transformation of the length x to radians, arc cos — or 
arc cos | ^ 


X a series of parallel cutting planes tangent to the Mach cone 

(At M = 1.0 these planes are perpendicular to the x axis.) 

Tlf angle in the xy plane forward between the y axis intercept 

of the cutting planes X on the xy plane, 

arc tan (-Vm®-! cos 0) 

^n(S) = Tchebichef polynomial, reference 6 

Vn(i) = = l^-l(g) Tchebichef polynomial, reference 6 


COMPUTIISG METHOD 


A s umm ary of the computing method is presented here, with the complete 
details given in the appendix. As shown in reference 1 (based on the 
theory of ref . 5 ) , the wave-drag equation may be written in coefficient 
form as 




n=i 


nA^^d0 


( 1 ) 
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wiiere A^ Is considered as having a length dimension and is defined as 


^n 


2 

1? 


1 


It 


S’ (x) sin nq) d 9 


( 2 ) 


Introduce the variable, g = ^ = cos q)j then ' equation (2) may be written, 

L 


An - 


n 


f\.( ) 

j ^ sin cp 


(3) 


where the limits are defined in figure 1: 

2 ^ ^ ^ 2 
- 1 < I < 1 

Jt > cp > 0 


Introduce the Tchebichef polynomial, reference 6: 


Vn(l) 


sin nq) 
sin cp 


Then equation (3) for the coefficients of the Fourier sine series may be 
written such that the integrand is a fvinctlon of one variable: 


An 



S'(|)Vn(l)d|, 


n=l,2,3. 




Because of this simplification the computation of these coefficients and 


the summation ^ --n 2 

n=i 


I 


can be performed by one continuous operation on 


a digital computing machine. Likewise, a reverse check computation can 
be performed by one machine operation. The wave-drag coefficients are 
computed from equation (l) by a simple integration as was done in 
reference 1 . 


There are definite advantages whicdi are characteristic of the new 
ccamputing method. It eliminates the intermediate steps of the previous 
method (ref. l) consisting of computing the slopes S’(x), plotting S’(x) 
as a function of cp, and then reading this slope curve. The new method 
works directly from the area-distribution curves of the model. An addi- 
tional .machine operation has been programmed (see the appendix) whicdi 
permits a chedk ccmiputation from the An coefficients back to the area 
curve. This is an all-inclusive check which may also be used to evaluate 
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•fctLe adequacy of tlie selected number of harmonics used in representing a 
particular cuirve. Supervision of the computation is cut to a minimum by 
the previously mentioned improvements. 

The machine time required to .make computations by the new .method is 
one-half that of the previoiis method (ref. l) . This comparison does not 
include the time lost using the previous .method due to data handling 
between steps . The time required to determine the area distribution of 
the models is not considered in this report and would not vary between 
the wave-drag computation methods . To compute the Fourier sine series 
solutions (25 harmonics) from pimched data cards representing one area 
curve ^ using a Magnetic Drum Calculator^ only 5 minutes are required. 

The check solution requires about 10 minutes . Thus on the assumption that 
five area curves are required per Mach number computation^ the machine 
time required for the Fourier solution using the Improved method would be 
less than I /2 hourj and the checking time, less than 1 hour. 

A complete derivation of the new computing method and the chec kin g 
procedure are given in the appendix which contains: 

1. Transformation to Tchebichef form for computing Fourier coeffi- 

cients 

2. General integration procedure 

3 . Tchebichef integration coefficients 

(a) Linear approximation 

(b) Qmdratic approximation 

i)- , Checking procedure 

5 . Construction and checking of tables (available on pmched data 

cards upon request) 

6 . Theory and Properties of Tchebichef Polynomials 

RESULTS AND DISCUSSION 


An evaluation of the new computing method will be discussed in three 
parts: known analytical solutions compared with computed values of the 

Fourier coefficients; previous solutions from reference 3 compared with 


new computed values 


of the drag parameter 


N 

n=x 


and available experi- 


mental values of drag-rise coefficients compared with computed values of 
wave-drag coefficients at zero lift. 
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Check of Method by Analytical Solutions 


The first known solution considered will be that for a Sears-Haack 
body. G?he shape of this .minimum-drag body for prescribed volume and length 
is defined in figure 2. The fineness ratio of 12.5 "the actual dimen- 
sions were arbitrarily selected. The theoretical equation for the zero- 
lift wave-drag coefficient (Cj) fig. 2) is seen to be independent of 
Mach number (for slender bodies or for Mach nimbers near l.O). Also note 
that the derivative of the area curve is exactly equal to Aasin 2cp where 
Ag defines the magnitude of the sin 2cp curve. The coefficient Ag is 
the only nonzero termj therefore equation (l) simplifies to 




t 



(5) 


where, in this case, the reference axe& (S^) is replaced by the body maxi- 
mum cross-sectional area for a body-alone test; that is. 



and 

^ nAa^ = 2As^ = ^ = 206. 521 sq in. 

n=x 

For the dimensions selected, d^ = 53.9198 inches. 


( 6 ) 


For comparison with this theoretical solution, the area curve in 
figure 2 was analyzed by the new method using Tcheblchef polynomials . 

The areas were computed for 200 equal increments of x, data cards were 
punched, and the A^ values were computed for 25 harmonics with the 
res^ilts as shown in table I. Even for the linear approximation (see 
appendix) between data points of S(|), the value of 2Az^ equal to 206. 502 
is in excellent agreement with the theoretical value of 206.521, and the 

2 5 


total summation 



equal to 206.731 shows an error only on the order 


n=i 


of 1/10 of 1 percent . For all harmonics other than the second the coef- 
ficients should be zero. The areas from the check solution agreed well 
with the original values for this case, with errors less than the possible 
errors in the original data (O.O5 percent of B(|) maximum). For a body 
that is closed at both ends the first coefficient of the Fourier sine 
series is always equal to zero, as in the solution presented in table I. 
For a body which does not close at the ends, the first term lAx^ represents 
the function of wave drag for a von KArmAn ogive on an infinite cylinder. 
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The cylinder cross-sectional area is equal to the base area of the body 
or is equal to the base area minus inlet area in the case of a ducted 
model. 

The second kno-wn solution considered will be that for a simplified 
area curve with a discontinuity of slope. This curve was selected in 
order to check the over-all computing procedure for all values of n £rom 
1 to 25 , because the previous Sears-Baack case primarily checked the compu- 
tation of the second harmonic. The area curve is shown by the small insert 
in figure 3 (a) and is defined: 

S(|) = 0 for -0.6 < I < 1.0 

and 

S(|) = -(1 + 0.6) for - 1.0 < I < 0.6 

The coefficients •A-n were computed by hand from the relationship 

-o .6 



-1.0 

These hand-computed values of An are compared with machine computations 
based on the area curve read at 200 equal intervals (table II) , Note that 
the machine solutions agreed very well with the hand computed values 
(within five decimal places) . These machine- computed values of An 
were used to compute the check solutions shown for both the area curve 
(fig. 3(a)) and the derivative of the area curve (fig. 3(b)). The area 
curve is enlarged in the region of the discontinuity of interest (| = - 0 . 6 ) 
in order to show the check solution differences (fig. 3(a)). At this dis- 
continuity the maximum error is 2.3 percent of S(|) .maximum^ otherwise 
the error is less than 1 percent. This type of check solution gives an 
obvious geometric plctiire of the accuracy of the computation for the number 
of harmonics used relative to the accxaracy with which the original area 
curve is computed. 

Although not required by the new computing method, the derivative of 
the entire area curve is shown as the solid line in figure 3(t>) . This 
slope curve is plotted against cp in the manner used in the previous 
method of computing wave drag (ref. 1). The circular check points were 
computed by the method of reference 3, but the same values of An (based 
on Tchebichef polynomials, fig. 3(s-)) were used. The overshoot of the 
check points near the discontinuities is due to the Gibbs phenomenon 
(ref. 8 ) and is not timly an indication of the accuracy of the Fourier 
coefficients , An- The interpretation of this type of check solution, rela- 
tive to errors in area, is not obviotis and is certainly more difficult 
than the check of figure 3 (a) . This difficulty was also apparent for some 
of the check cases presented in reference 3 - 
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Clieck of Method by Prior Solutions 


Although more involved, the previous computing method has been fairly 
successful in evaluating wave drag; therefore, a direct comparison between 


the two methods of the solutions of the drag parameter 



was made. 


Prom reference 3, two area curves were selected which had similar 


total summations ^ iiA^ , one curve having a discontinuity in slope of 
n=i If 

the area curve (figs. 4(a) and 4(b)). The drag parameters ^ nAn^ 


n=i 


computed by the two methods are compared in table IH. The agreement 
Indicates that the Tchebichef solutions are comparable within 2 to 3 per- 
cent even for a cinrve with a discontinuity of slope. Due to the reduced 
number of operations Involved, the new method has less chance for error . 


The check solutions, of these two area curves shown in figure 4, 
were computed by the new method, using the Tchebichef values of Aj^ 
(table IZE) . As was noted in reference 3 , the check solutions indicate 
a fairly satisfactory fit of the area ciirve in the case of figure 4(a), 
and an incomplete solution for N = 25 in the case of figure 4(b) . How- 
ever, in contrast with the check solutions of reference 3 > these check 
solutions directly show the errors in area. 


Check of Method by Experimental Results 


Available geometric and experimental data for three models were 
utilized to check the new method: a triangular-wing interceptor-type 

airplane model, a swept-wlng interceptor-type model, and a free-flight 
.model having tall surfaces but no wing (this model was also considered 
with a large scoop- inlet and duct system) . 

Triangular-wing interceptor-type airplane model - Model 1 .- The 
dimensions and pertinent geometric data of the model are presented in 
figure 5 (data from ref. 9 ) . The experimental data, also from reference 9^ 
are shown in figure 6 with the computed results and model cross-sectional 
area distribution. The experimental drag coefficients are identified as 
minimum-drag coefficients, because the lift varied from approximately 
zero lift at low subsonic and hl^ supersonic speeds to a peak value of 
C-Q = 0.05 at M = 0 . 97 . The friction-drag coefficients were approximated 
by the method used in reference 10. The empirical correction factor of 
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1.15 siiggested (ref. 10 ) for the body portion was not applied since good 
agreement was obtained with the subsonic drag coefficients without such 
Correction . 

The agreement between the experimental and theoretical total-drag- 
coefflclent curves (fig. 6 ) is considered to be good, with errors less 
than 10 per cent \of the experimental minimum-drag coefficients. The error 
in percentage of the wave-drag-co'efficient increment, assuming the 
ffiction-drag variation is correct, would be increased to almost I 6 per- 
cent. An assumption of no variation in £rlction-drag coefficient with 
Mach number would indicate essentially no error at all supersonic speeds . 

As was discussed in reference 3 , the theory does not apply when the leading 
edge of a wing with a roxmded leading edge becomes sonic. For this con- 
figuration the leading edges of the horizontal stabilizers are supersonic 
for Mach numbers greater than approximately 1.4 and the wing leading edge 
is supersonic for Mach numbers greater than approximately 1.6. Inaamuch 
as the wing is thin and the horizontal tall is a small part of the total 
configuration, the computations at M = 1.8 are considered to be a fair 
approximation in spite of the violation of the assumptions of the basic 
theory. The leading edge of the veirtical surface of the tail is not sonic 
until M = 2.0. 

As a secondary evalmtion of the computations for this triangular- 
wing model, area curves with check solutions (IT = 25 ) are presented in 
figure 7 for Mach numbers of 1.0, 1.2, and 1.8. These cxjrves illustrate 
the degree of eiror in the area curves which are apparently tolerable, 
in view of the good computed drag- coefficient results shown in figure 6. 
This suggests that the area curves are effectively smoothed by boundary- 
layer or separation effects. For the M = 1.0 curve (fig. > ccmputing 

additional terms of the series would probably Improve the agreement between 
the check curve and the original curve, but would Increase the disagreement 
with the drag-rise coefficients at M = 1.0. The check solutions for the 
M = 1.2 curves (figs. 7(b), 7(c), and 7(<3.)) considered to be satis- 
factory because only slight radii changes to the original body would be 
required to give areas which would .match the check curve. Experimental 
results have Indicated that small changes in area of this order or less 
would have little effect on the wave drag. The check solutions for the 
M = 1.8 curves (figs. 7(e) en<i 7(g)) ere considered on the basis of the 
area curves alone to be questionable, particularly for the 270 ° cut 
(fig. 7 (g)). 

In this case, for the M = 1.8 ccmputatlons , additional haimonics 
over the 25 used are evidently required to Improve the agreement between 
the area curves . Increasing the harmonics to ^9 gave agreement between 
the area cinrves and the check solutions (generally within the width of 
the line), with the maximum error ( 270 ° cut) being less than 2 - 1/2 percent 
of the local area curve values. The computed wave-drag coefficients for 
IT = 49 were 8 percent greater than the experimental drag-rise coefficients 
at M = 1.8 (fig. 6 ) . Increasing the number of harmonics above i^9 would 
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further improve the agreement with the area curve, but would increase the 
disagreement with the experimental drag resiilts . This tends to support 
the previously mentioned fact that at M = 1.8, the supersonic leading 
edges of the blunt airfoil sections are in violation of the theoretical 
limitations of smoothness or slenderness. 

Swept-wing interceptor-type airplane model - Model 2 .- The dimensions 
and pertinent geometric data of this model are presented in figure 8. The 
dimensions shown were taken from layout drawings furnished by the manu- 
facturer for this particular version. The layout drawings also contained 
details as to fillets, etc., from which comprehensive area distributions 
were computed. Flight-test resiolts were obtained from reference 11 for 
a model of a similar low- tail version of the airplane. The area distri- 
butions for Model 2 and the model of reference 11 are compared in figure 9> 
The essential difference between the models is that the low-tall version 
of reference 11 had thinner tall surfaces. In addition the fuselage lines 
on the model of reference 11 were faired to a point ahead of the proposed 
nose- inlet location and the model had no duct. 

The two area curves shown in figure 9 were analyzed by the new com- 
puting method (N = 25 ) which gave wave-drag coefficients for M = 1.0. 

The computed resilts are compared in figure 10 with experimental drag-rise 
coefficients from reference 11. The only exact comparison which can be 
made is for the experiment and computation for the model of reference 11 
for M = 1.0, The perfect agreement is fortuitous, particularly at a 
Mach number of 1.0, but still this is a favorable indication for the new 
computing method. The difference between the computed resilts for the 
two models (M = 1.0) is only sli^tly over 10 percent, based on the experi- 
mental results. 

The computations at the two higher Mach numbers had to be made for 
the airplane for which detailed area distributions were available; there- 
fore, a direct comparison between the theory and the experiment cannot 
be made. However, the differences at these higher Mach numbers should 
be of the order of the differences at M = 1.0 (or less), thus reasonable 
agreement is indicated. No correction for friction drag was made in this 
case; however, the prior example showed that the variation in frlctlon- 
drag coefficient with Mach number would be slight. Computations were not 
attempted above M = 1.4 becaiise sonic leading edges of the fairly thick 
wings with roimded leading edges would produce abrupt changes in area 
distribution which definitely exceed the theoretical limitations . For 
this type of wing the 25 "term solution does ^ot closely fit the original 
area curve in the region of the discontinuity, as shown in reference 3 
and figure 4(b) . 

Scoop-lnlet-duct model - Model 3 .- The results of the experimental 
investigation of this model have not been published; however, the test 
procedures were identical with those of reference 12. Details of the 
models and a brief discussion of the experimental data accuracy are 
Included, so that comparisons with wave-drag coefficients computed by the 
new method may be correctly evaliaated. The dimensions of the basic body 



NACA EM A55J28 


11 


and the ducted version are given in figure 11 and tables IV and V . The 
tests of the scoop-duct version were tiiade at duct tnass-flow ratios of 
approxitnately 0.6 and 0.9 using two duct nozzles with different throat 
sizes . The accuracy of the instrumentation and the accuracy of the final 
data were similar to the values of reference 12 . For the scoop model , 
particularly for the test with a mass-flow ratio of 0.6, the errors at 
subsonic speeds were increased, due to fluctuation of duct pressirres , to 
a possible Tnaximum error in drag coefficients of ±0.02. 

The experimental results for Model 3 s.re presented in figure 12. 

The internal drag coefficients (computed by the method used in ref. 12) 
were subtracted from the total drag coefficients at each data point to 
obtain the external drag coefficients . External drag-rise coefficients 
were measured relative to the experimental data at a Mach number of 0.8 
(fig. 12(b)). The small differences between the experimental drag-rise 
coefficients at the two different mass-flow ratios are probably not sig- 
nificant in view of the possible errors in these data at subsonic speeds . 
The external drag-rise coefficients of the basic model are of course more 
accurate than for the ducted models, because the basic data were more 
acciirate and no correction for internal drag was required. 

Computed wave-drag coefficients are also shown in figure 12(b) . The 
computation for the basic model was stral^tforward, checks were satis- 
factory, «-nd the agreement with experimentation can be considered as 
another check of the new method of computation (differences k 10 percent 
or less). In the computations for a mass-flow ratio of 1.0 it was assumed 
that a streamtube eqioal to the duct-inlet area passed through the duct 
relatively undisturbed without contributing to the external drag coeffi- 
cients (fig. 13 ) . The duct- inlet area was removed as a constant value 
from fuselage station 8I.5 (duct lnlet)to station 13 ^ (duct exit). The 
area difference between the duct inlet and the exit was added as a constant 
value to the ordinates of the portion of the model aft of the duct exit , 
on the assumption that the increased diameter of the streamtube is main- 
tained and affects somewhat the external wave drag. Another simplified 
procedure for determining the equivalent area distribution was tried which, 
in this case, gave essentially the same results. The area removed was 
determined by taking a straight-line variation from the inlet area to the 
exit area. As in the prior method, this procedure removed the steps (at 
the inlet and exit) in the original area curve. 

The computed wave-drag coefficients for (mi/mo) = 1.0 were compared 
with experimental drag-rise coefficients for (mi/mo) ~ 0.9^ inasmuch as 
the experimental variation of external drag-rise coefficients with mass- 
flow ratio is slight at transonic speeds (e.g., see data of ref. I3, con- 
verted to drag-rise coefficients). Agreement between experimental and 
computed values is good with max^Tnum differences of the order of 6 percent 
of the experimental values. This example is not considered as a firm 
indication of the accuracy of the computing method because the experimental 
data co\ld be in error by 10 percent. Within the evident accuracy of the 
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theory and experiment, the check solution for N = 25 (fig. 13) Indicates 
an adequate representation of the area curve. Additional terms would be 
required for an exact fit. 

The computations for a mass-flow ratio of 0.6 used the latter 
described procedure for removing area. Removing an area varying linearly 
from 0.6 of the inlet area to 0.6 of the exit area left sizable steps in 
the area curve (solid line, fig. 1^). These large steps in the area curve 
result in infinite slopes and a configuration which is outside the Intended 
limits of the theory. Some arbitrary fairings could be inserted at the 
duct inlet based upon estimated location of the shock wave ahead of the 
duct, but even this fairing would be rather abrupt and would not remove 
the discontinuity at the exit. Additional theory and experimentation are 
required to know how to correctly handle this type of discontinuity. In 
spite of these limitations, a 25 -term computation was made to clarify the 
difficulties and to Illustrate why the computation is an inadmissible one. 
The computed wave-drag coefficients for a mass-flow ratio of 0.6 predict 
a large increase in drag-rise coefficient with a decrease in mass-flow 
ratio from 1.0 to 0.6. However, as in reference I 3 , the experiments indi- 
cate little variation in drag-rise coefficient with mass-flow ratio. 

The check solutions for N = 10 and H = 25 (fig. l^(a)) indicate that 
the 25 “term solution is not adequate to fit the curve and at least double 
the number of terms would be required to even approach a close fit. 
Increasing the number of terms above 25 would of course increase the dis- 
agreement with measured drag coefficients, and this case is clearly an 
inadmissible computation.^ At present it is recommended, for ducted con- 
figurations, to make the computations of wave drag only for a mass -flow 
ratio of 1.0, and estimate the possible variations with mass-flow ratio 
from experimental results. 


CONCLUDIHG RIMARKS 


The use of Tchebichef polynomials simplifies the computation of zero- 
lift wave drag and provides a direct check of how well the number of 
harmonics used represents the area- distribution curve. For relatively 
smooth distributions of area, about 25 harmonics are usually sufficient 
to represent the area curve and to compute wave drag in reasonable agree- 
ment with experimental data. 

Even for ccmpl^ airplane configurations with sll^t discontinuities 
in slope of the area distribution, 25 harmonics gave reasonable agreement 

^It was of Interest to note that a very arbitrary selection of 
18 terms gave drag results in agreement with the experimental data, but 
the check solution for I 8 terms (fig. l^(b)) resulted in an extreme fair- 
ing of the area curve. 
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with the measured drag rise although the calculated values were for a 
smoothed area distribution. This suggested that small discontinuities 
might actually be smoothed by boundary-layer or separation effects. 

On the other hand, for configurations with large discontinuities in 
slope (e.g., produced by a duct with a mass-flow ratio other than l.O), 
the computed "smoothed" area distribution and wave drag varied .markedly 
with the number of harmonics. This is to be expected because if an infi- 
nite number, of harmonics were used, an infinite drag would be predicted 
when there is a step in the area curve. 

In contrast a smooth body may require very few harmonics (as low as 
1 or 2) to give an exact fit of the area curve and an acciirate wave-drag 
prediction. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Oct. 28, 1955 
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APPENDIX A 

TCHEBICKEP FORM OF HARMONIC AMLYSIS 
Transformation to Tcheblclief Form for 
Computing Foirrier Coefficients 


In reference 5 the strearawise variation of cross-sectional area 
S(x) is expanded in a Fourier sine series in cp: 


M 

dx 


00 

= ^ Ansin n£p, 
n=i 



<1 


(Al) 


where the Glauert angle, cp, is given by 


cos cp 



0 < cp < :ir 


(A2) 


The basic problem then is to devise a practical procedure for computing 
the Fourier coefficients, An. To begin, they are given by the classical 
foimula of harmonic analyses ; 


An = ^ y ~ sin ncp dcp (A3) 

To transform to the Tchebichef form, introduce the dimensionless coordi- 
nate, 

I = = COB cp, -1 < I < 1 (Ai^) 

1/2 - - 


n=0,l,2. 


Introducing the Tchebicihef polynomials, reference 6, 

1 

Tn(l) = cos ncp 

^n(l) Wl), 

one may write eq^uation (A 3 ) as 


n=l,2,3, . . 


(A5) 


(A6) 


-1 


nn NFTnENITTAT. 
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(a summary of the theory and properties of the Tcheblchef polynomials is 
given at the end of this appendix.) Dae to the fact that 

# equation (A^tO), equation (a 6 ) can be simplified to 


An = 


k_ 



^n(S) 


+ (-l)^S(-l) 




(AT) 


Since Vi(|) =1, from, the definition, equation (A 5 ) , the first coeffi- 
cient, Ai, can be computed Imnedlately; 


Ai 



S(l) - S(-l) 


(A 8 ) 


Henceforth A^ will be considered only for the positive Integral values 
of n = 2 , 3 , 5, . . . . 

In the next section general numerical Integration formulas will be 
developed, and in the following section they will be applied to 
equation (AT) . 


General Integration Procedure 


The general problem of the numerical evaluation of integrals of the 

type 


f 


f(l)g(l)a| 


has been dlsciissed by Sheldon in reference l4, and the particular case 
in which g(|) is a trigonometric function has been treated by Filon in 
reference I 5 . Here the problem is generalized to include any derivative 
of f(|) rather than just the function Itself. The so-called "strip 
method" is used, in which the range of integration is broken into H 
strips of equal length, and f(|) is approximated in each strip by a poly- 
nomial . Thus 


H-i 


®-k+i 


/ u T n 

f(E)g(S)<ig =2, J fCDsCD'il 

k=o 


(A9) 


where 


= a, 




a|j- = b are the boundary points of 


the If strips. The integral over each strip is evaluated by approximating 
f(l) t)y S’ polynomial of degree m and sxibdividing the strip into m 
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intervals of length A, this interval length being constant throu^out 
the entire range (a,b) . Thus, the lattice points within the kfch strip 
are aj^, ajj. + A, . . aj^ + .nA = Since the case of most interest 

is that in which f(|) represents experimental data, it is assumed that 
f(|) is not known continuously but only at these lattice points. Hence 
it is desired to obtain an integration formula involving only these 
points; 

'^-k+x 

(AlO) 

/ t 


J f(&)g(S)d.g =2^ a(m,j,k)f(a]s. + JA) 

J=o 


where the coefficients a(m,j,k) are to be determined as functions of 
g(|). Substituting eq-uation (AlO) in (A9) gives 


H-i m 


J f(§)s(S)<iS = ^ a(m,j,k)f(ajj. + 




k=o j=o 


(All) 


If all the lattice points are numbered consecutively over the entire 
range, a = = b, then 


5km+J - ®-k + 


and rearrangement of equation (All) gives 



Hm 

f(l)g(g)<a6 = 

r=o 



li(m,r)f(|j.) 


(A12) 


where the t-L*s are just the a’s suitably renumbered, except that the 
IJ. corresponding to each strip boundary point contains two a’s, one 
being the last a for the left-hand strip and the other being the first 
a for the right-hand strip, thus; 


p,(m,0) = a(m,0,0) 
p.(m,km) 


a(m,0,k) + o(m,m,k-l) 
k=l,2,3, - . lT-1 


n(m,km + j) 

M.(m,Nm) = a(m,m,ir-l) 


a(m,j,k), k=0,l, . . ., N-1 

j— 1,2, . . •, m— 1 


> 


(A13) 
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TtLUS equation (A12) is the desired integration formula, the coefficients 
IJ. being given by equation (AL3), while the coefficients a £tre implicitly . 
defined by equation (AlO) . It remains to determine them explicitly. This 
is done by requiring that equation (AlO) be exact whenever f(|) is a 
polynomial of degree m. (it must be remembered that all numerical inte- 
gration f03cmulas are- approximations when the integrands are arbitrary.) 

This is equivalent to requiring that equation (ALO) be exact for any set 
of m -t- 1 linearly independent polynomials of degree m. The most conven- 
ient choice is the set of Lagranglan interpolation polynomials L(m,J,kj|) 
defined by (ref . l6) : 


L(m,j,kja]5. + lA) 


0 if X ^ i 

1 if r = J 


(ALlt) 


An explicit representation is 


L(m,j,kj|) 


m i - (a]j: lA) 

i?o »k + “ (^■k + 


j‘.(m - j)lA’^ 


m 

n 

1=0 


“ ^km+i^ 


i/j 


(AL5) 


Substituting L(m,J,kj|) for f(|) in equation (AlO) gives, by virtue of 
equation (Al^t) the explicit formula for the a: 


/->®-k+x 

a(m,j,k) = / L(m,j,k;|)g(|)d| 


(Al6) 


not 


More generally, if it is desired to evaluate an integral involving 
f(|) but any one of its derivatives, 

/ 


r'‘ fm. 


g(l)<ii 


then formulas (AlO) throu^ (Al6) remain valid, except that in (A16) L 
must be replaced by the corresponding derivative 


d^L(m,j,k;g) 

d|^ 

while f is replaced by the derivative in the integrands but not in the 
sums . 
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Clearly tlie effectiveness of this entire procedure hinges upon the 
dlfflctilty of calculating the a from equation (Al6) . To estimate this 
difficulty It Is worthwhile exhibiting some of the L’s expillcltly; the 
second form of equation (AI 5 ) gives, for m = 1, 2, 3^ respectively: 


L(l,0,kjg) = - | (g -A+i)' 
L(l,l,k;|) =s (S - gfc) 


L(2,0,kj|) 2^a (S “ “ ^ak+a^ 

(1 ” lak^^^ “ ^ak+a^ ' 


L(2,2,k;g) = 1 ^ (g - gak)(S - ^ak+J J 

L(3^0^1t>g) - " ^3k+x^^^ “ ^ak+a^^^ “ ^sk+a^ 

L(3^1^J5:j|) = 2 ^^ (g - ggjj)(g - ” ^3k+3^ 

L(3^2,kj|) = (g ” gak)(^ " ^sk+i)(5 " lsk+ 3 ) 

L(3,3,kjg) = (I - S3k)(S - WHS - gskn^) . 


(A 17 ) 


(A18) 


(AI9) 


It Is a stral^tforward, thou^ tedious, task to extend such a table 
to as high a degree, m, as Is desired. The calculation of the Integral 
In equation (AI 6 ) then depends on the nature of the function g(g). This 
calculation can be performed analytically whenever g(g) Is a polynomial, 
trigonometric, hyperbolic, exponential, or logarithmic function. For 
more complicated functions, g(g). It may be necessary to resort to numeri- 
cal Integration at this stage, which almost defeats the purpose of the 
entire procedure! However, these Integrations (eq. (AI 6 )) are performed 
only once, and then equations (A12) and (AI 3 ) are applicable to arbitrary 
functions f(g). In the next section, then, these results will be applied 
to the Integrals of equation (A7) . 
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Tcheblclief lategration Coefficients 


Writing equation (Aj) in the form, of the preceding section^ 


Aji - 





Vn(l)i| 


(A2D) 


■with 


Sn(l) = S(g) + 

the desired resialt is, analogous to equation (A12) , 

M 

A-n = M-nr^n(Sr^ 


(A21) 


(A22) 


where M = mlT, the subscript m has been dropped for bre-vlty, while the 
subscript n has been added for cleirity. The p*s are still given by 
equations (AI 3 ) , but (AI 6 ) now becomes 


a(m,j,k) = / 
"4k 




’^n(E) 


dL(m,j,k;g) 

d| 




(A23) 


These ■will now be integrated explici'ty for the linear (m = 1) and -quad- 
ratic (m = 2) cases . 


Linear approximation .- In the linear case each strip consists of a 
single interval of width A, m = 1, a^ = = kA, M = N, NA = 1, and the 

derl^vatives of the Lagranglan inteirpolation polynomials are, from 
equation (AI 7 ) 


dL(l,0,k;S) ^ 
Substitution in equation (A23) 

a(l,0,k) 

a(l,l,h) 


^ , dL(l,l,k; I) ^ 1 
A ' d| A 

gi^ves, for j = 0 , 1 , respectively. 








^k+i 


> 


i / " Vt)ii 

^k 


(A24) 
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- jJJl 

Wow, by equation (A¥f) ^ ^ therefore equation (A2^) can be 

integ3:*ated directly to give 


a(l,0,k) 




a(l,l,k) 




Substituting these in equation (AI 3 ) and changing the notation to agree 
with equation (A22) gives 


4 1 

^0 “itnA 

Tii(O) - Tii(A) 




M'nr = - — r 
jtnA 

[-TnUr-J + 2T^(6r) * Tn(lr+x)j 

/ r=l,2. 


“ 

Tn(l) - Tn(l " 

A) 






* j 


N-1 


> (A25) 


In the linear case, then, the ccaapCLete solution is given by equa- 
tions (A21) , (A22) , and (A25) • 

Quadratic approx^wlAtion .- In the quadratic case, each strip consists 
of two Intervals, each of width A, m = 2, M = 2N, a^^ = = 21sA, 2HA = 1, 

and the derivatives of the Lagrangian ihterpolatlon polyncanlals are, from 
equation (AI 8 ), 


dL(2,0,k,-g) 

dg 

dL(2,l,k; g) 
dg 

dL(2,2,k;g) 


2g - (4k + 3)A 
2A^ 

(4k + 2)A - 2g 
A2 

2g - (4k + 1)A 

2A^ 
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Su'bs'fci’fcu'bioii in equation (A23) gives, for j - 0,1,2, respectively. 


■2 21 - (4k + 3)A v^(|.)d| 
2A2 


a(2,0,k) = j 

^2k 

0(2,2,.) ./^^-St^i^±iAv„(t)dt 




(A26) 


Now, from, equation (a42) , 

2|Vn(|) = Vn+i(|) + 

and, as before, i If these relations are substituted in equa- 

tion (A26) , the integrations can be performed lumediately to give 

a(2,0,k) = w(n + l,k) + W(r - l,k) - (4k + 3)W(r,k)Aj 
a(2,l,k) = ~ (4k + 2)w(n,k)A - W(n + l,k) - W(n - l,k) j 

a(2,2,k) = ^ W(n + l,k) + W(n - l,k) - (4k + l)W(n,k)Aj J 


where 


W(a,k) = 


(A2T) 


(A28) 


For this case equation (AI 3 ) becomes , in the notation of equation (A22) : 

l%,o = o'<(2,0,0) 

L, 2 T “ '5 “>(2^0^3r) + a(2,2,r - l) 




to,2N ~ ^ a(2,2,N - 1) 


(A29) 


Tn the quadratic case, then, the complete solution is given by 
equations (A21) , (A22) , and (A27) to (A29) . 
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It Is clear ftom eq^uations (AI 5 ) and (A23) tliat in order to use 
higher degree approximating polynomials, it is necessary to evaluate inte- 
grals of the type 


f I'^VnCDdl 

1 km 

for successive values of j = 0, 1, 2, . - m - 1. While no explicit 
formulas will be given, the section "^eory and Properties of Tchebichef 
Polynomials” does Include some recurrence relations which can be used in 
an obvious way to perform such integrations systematically for successive 
values of j . In every case the results, a and ti. Involve only the one 
kind of Tchebichef polynomials Tn(l)^ ^n±i(&)> etc. 

Numerical tables of the integration coefficients lijQQ. have been 
constructed, by the procedures discussed below, for the linear and quad- 
ratic cases. 


Checking Procedure 


In any lengthy numerical calculation it is desirable to have a check- 
ing procedure, even if it merely consists of repeating the calculation. 

In the present case it is possible to do much better than this. Once the 
Fourier coefficients, Aq, have been calculated, they can be used in a 
reverse sense to calculate the area distribution, S(|), which can be com- 
pared with the original data as a check. This is a quite effective check, 
since it detects any of the following; eirors in the theory Itself, errors 
in the tables of integration coefficients, errors in reading graphical 
data, errors due to using too few terms of the Fourier series (eq. (A2)), 
round-off errors, and errors made by the computing machine or human opera- 
tor. Clearly, such an independent checking procedure is preferable to a 
mere repetition of the direct ceilculation, which would only detect machine 
or operator errors. 


Such a check is easily derived from equation (A2) . Transforming to 
the variable | by means of equations (Ak) and (A 5 ) gives 


i'll 

n=i 


From equation (a45) , 




f 





Vn+i( £) 
n + 1 


Vn-i(l)l 
n - 1 ' ' 


(A 30 ) 


n=2,3. 
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while, since Vi(|) = 1, 


ViCD-TiTF 



+ ~ arc sin | 


Substituting these .in equation (A 30 ) gives 


M = 1 A. 

d| 4 d| 



CO 


+ arc sin | )Ax + ^ 


n=2 


'Vn+xU) 
. n + 1 


Vn-i(g) 
n - 1 


and S(|) can be obtained immediately. 


S(l) = C + -^ + arc sin I^Aj. + i ^ 

n=a 


where C is a constant of integration. Setting | = 1 gives 

S(l) = c + ^ Ai 

Eliminating C in equation (A 32 ) gives 


SU) = S(l) + 1 (1*71^ - arc cos 


00 

/ ^ L ^ 1 


It is also convenient to set | = -I in equation (A 32 ) , giving 


• ' 

« 

(A3I) 

x(g)~ 

- 1 _ 

(A32) 


Yn-x(g) 
n - 1 

(A33) 


S(-l) = C - ^ Ai 



2h 
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Using this to eliminate C in equation (A 32 ) and then substituting -§ 
for § gives 


S(-|) = S(-l) - 1 



- arc cos 


1 

4 



(-1)% 


Vn-n(^) 
n + 1 


+ 

Vn-iU) 

n - 1 


(A34) 


by virtue of equation (a 40) . Equation (A33) is used as a check for posi- 
tive values of | and (A 34 ) for negative v^ues . Ehey can be -written 
more briefly as 


00 

S(±|) = S(±l) + I (±l)“-Cn(i)A^ 0 < s 1 


n=i 


•where 


(A35) 


Old) =|(en/i^ 
Cn(g) 


„ arc cos | 


Vn+x(g) Vn-x(g)' 


n + 1 


n - 1 




(A36) 


and in (A 35 ) either the upper or lo-wer sign is to be used throughout. It 
should be emphasized that equations (A35) ^<3. (A 36 ) constitute a check 
regardless of -the choice of any particular method of computing the Fourier 
coefficients , An • 


Construction and Checking of Tables 


The usefulness of the present theory lies in the fact that it makes 
possible the construction of numerical tables of -the Integration coeffi- 
cients, equations (A25) and (A29), and of the check coefficients, 

Cn(|)^ equation (A 36 ) . Once such tables have been prepared, the calcula- 
tion of the Fotirier coefficients, Aq, fSrcraa equation (A22) and the checking 
by means of equation (A 35 ) is a straightforward, systematic, arithmetical 
procedure, whether or not automatic computing machinery is a-vallable. Such 
tables of coefficients have been ccaaputed (and are a-vallable on punched 
data cards upon request) for Integral -values of n from 1 to 98 , and for 
I in -the range of 0 to 1 at an inter-val of 0.01. All these calculations 
were performed on an IBM type 605 Electronic CalcTxLatlng Punch. Fixed- 
decimal, 13 -digit arithmetic was used, and the final -values of the 
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coefficients and Cn(|) were rounded off to 10 digits to be in a foim 

suitable for use in magnetic-drum calculators. The sequence of calcula- 
tions was as follows (see next section for proofs); 


a. To(g) = 1, Ti(|) = I 

~ 2|Tjq(|) — n — l, 2 , 3 j ***^99 

c. Vo(|) = 0, ViU) = 1 

< 3 - Vn+i(§) = 2|Vii(|) - n = 1 , 2 , 3,-.., 99 

e. Qa(S) was then computed lislng eguatlon (A36) with n = 1 , 2 , 3 , • ■ -,98, 

I = 0.00, 0.01, 0.02, .. .,0.99, 1.00 

Tn(l) - Tu^(£ - 0.01) was computed for n = 0,1,2, ... ,99, I = 0.01,, 
0.02, . . .,1.00 

g. P]xr(n = 2,3, •••,98, r = 0,1,..., 100) was computed for the linear case 

from equations (A 25 ) • 

h, W(n,k) was computed for n = 1,2, ...,99, h = 0,1,..., ^9 from equa- 

tion (A 28 ) 

1. (x(2,j,k) was computed for J = 0,1,2, k = 0,1,2, ...,49 from equa- 

tion (A 27 ) 

j- fJ-nr(°- = 2 , 3 , ••.,98 j r = 0 , 1 ,..., lOO) was computed for the quadratic 
case from equation (A 29 ) 


As a preliminary check, the tables of Tn(|) and Vq(|) were checked 
using equations (a43): 


Tn+i(l) = |Tn(l) - (1 - l^)Vn(l) 
Vn+i(l) = |Vi,(S) +Tn(|) 


n= 0 ,l, . . .,98 


The tables of were checked as follows; For the linear case, the 

integration formula (A22) shoiild be exact if S(|) is linear. That is, 
equations (A20) and (A22) should give the same resilLts for the Fourier 
coefficients. An, whenever S(|) is linear. This check was applied for 
various linear functions S(|). For the quadratic case, equations (A20) 
and (A22) should give the same result for the Fourier coefficients. An, 
whenever S(|) is linear or quadratic. This check was applied for various 
linear and quadratic functions S(|). 


The tables of check coefficients, Cn(|), were checked by construct- 
ing analytic functions S(|) having oifLy a finite number of nonvanishing 
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Pourler coef f icients , A^, a.i i of wiiicli were calculated analytically and 
then substituted in the check equation (A 35 ) . 

The purpose of these checks is to ensure that all the tables are 
free of numerical errors. There still remains the question of two other 
sources of error; the approximation of S(|) by a polynomial, and the 
use of only a finite number of terms Ajj^. These two effects are discussed 
in the main body of the present report, jerticularly in conaectlon with 
the Sears-Baack body. In general, such questions can only be examined 
empirically. If the check procedure indicates the presence of errors, 
three things can be done. More coefficients A^ can be taken, a finer 
interval A can be used, or higher degree polynomial approximation can 
be used. Any of these courses of action requires extending the basic 
tables. In many practical cases the data function S(g) is originally 
obtained by fairing a curve through a reasonable number of points . Obvi- 
ously, the use of a finer reading Interval A is incapable of yielding 
greater accuracy in such cases. Use of hi^er degree polynomials is some- 
times justified, birfc this also can give illusory improvements. For 
example, a curve with sharp changes in slope is better approximated by 
short straight-line segments than by longer arcs of higher degrefe poly- 
nomials. Finally, the use of more terms in the Fourier series can be 
reduced to an absurdity once the wave lengths involved became shorter than 
the reading interval A. 


Theory and Properties of Tchebichef Polynomials 


The Tchebichef polynomials are defined in reference 6 as 
Tn(S) = cos ncp 
Un(S) = [sin(n + l)cp]/sln 9 


(A3T) 


where | = cos 9. Since Un-i(S) is closely related to the Fourier coef- 
ficient, An, it seems convenient to introduce the notation 


Vn(l) = Un-x(l) = (A 38 ) 

Bln 9 

Since the principal properties of these polynomials. Including the fact 
that they are polynomials in are treated extensively in reference 6, 
the present section will list and prove only those that have been used 
in this report. 


First, the following are obvious from the definitions; 


To(i) 


= 1 , 


(A39) 


^x(S) = I, VoU) = 0, Vi(|) = 1 
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Next consider tlie effect of replacing | by so that 9 is replaced 

by It - cp, that is, cos(jt - 9 ) = -cos 9 : 


Tn("t) = cos n(ic - 9 ) = cos n% cos n 9 


y _ sin n(ic - 9 ) _ -cos n:t sin n 9 

^ sin(:t - 9 ) -cos sin 9 


Thus 


Ta(-S) = (-l)''Tn(|), 

for all integral values of n. 


(a4o) 


Next, setting g =1 gives 9 = 0, T^ = 1, while is the indeter- 

minate form 

llm 

0 sin 9 


9 

Applying L’ Hospital’ s rule gives 

Tn(l) = 1, Vn(l) = n 


(Al^l) 


Recurrence relations can be obtained directly from the definitions 
using the trigonometric addition formulas 

sin(a ± b) = sin a cos b ± cos a sin b 

cos (a + b) = COB a cos b ± sin a sin b 

By use of these, the following are easily proved: 

Tn+i(l) = 2 £Tn(|) - 

Vn+i(5) = 2|Yn(S) - Vn_i(g)^ 

and 

Tn+iU) = ETnd) - (l - . 

Vn+id) = |Yn(g) + Tnd) 


(a42) 


(a43) 


Numerous differential and integral identities are obtainable, analogous 
to the trigonometrical formulas 

d sin Tjr , d cos i]f 

~ dt = — d T = ^ 
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For example, differentiating the definition of £ gives 


— = -sin cp = 
dcp 

Thus, differentiating Tn(£), equation (A 37 ), gives 

<3-T]2 _ dTj^/dtp _ -n sin ncp 
d| d|/d<p -sin cp 


that is 


aind) 


= nVn(l) 


(a44) 


Next consider 


d^ 




Vn+id) Vn-i(l)' 


n + 1 


n - 1 


< 

d 

► ‘ 
1 fslnfn + l)cp sln(n - l)cp"l 

- 


2 \_ n + 1 n-1 J 

, 4 


= i j^cos(n + l)cp - cos(n - l)cpj || 

= y = VnBin cp 

-Bin cp n T 


that is 


A 

d-E 




Vn+iU) Va-iU)' 


n + 1 n - 1 
Finally, to evaluate integrals of the type (see above) 


= Vj,(|)n/iTF 


(a45) 


I(j,n) s J E^VnCDdi (a46) 

consider first the case J = 0. This can be integrated immediately by 
virtue of equation (aMi-) to give 

l(0,n) =iTn(£) (a4T) ’ 
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Next, If j > 1, use the recuirrence relation (a 42) to write 


I(j,n) = = Ij/" 


Vn+x(l) +Vi,.3.(g) 


< 3-1 


that is 




1 

2 


I(j - l,n - 1) 


+ I(J - l,n + 1) 


(a48) 


Thiis, from equation (a 47), l(0,n) can be evaluated for n = 1,2, . . . . 
Then from equation (a 48), l(j,n) can be obtained by recurrence, us ing 
the one additional equation 


since Vo(i) = 0. 


I(j,0) = 0 


(A49) 


By means of these equations the integrals l(^,n) can be evaluated 
systematically, either analytically or numerically. The general integra- 
tion procedtn-e presented above can then be applied to extend the Tchebichef 
integration coefficients to any desired degree of polynomial approximation. 
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TABLE I.- EAE?MOITIC SOLUTION FOR A FINENESS RATIO 12.5 
SEARS-HAACK BODY BY TCHEBICHEF POLYNOMIALS 


n 

Linear 

approximation^ 

An 

Lrfn^r 

approximation! 

N 

I 

n=x 

Quadratic 
approximation ! 
N 

z 

n=x 

1 

0.00000 

0.000 

0.000 

2 

- 10.16125 

206.502 

206.504 

3 

.00000 

206.502 

206.504 

k 

.00091 

206.502 

206.504 

5 

.00000 

206.502 

206.504 

6 

-.00091 

206.502 

206.504 

7 

.00000 

206.502 

206.504 

8 

-.00821 

206.502 

206.504 

9 

.00000 

206.502 

206.504 

10 

- .01982 

206.506 

206.505 

11 

.00000 

206.506 

206.505 

12 

-.02552 

206.514 

206.506 

13 

.00000 

206.514 

206.506 

14 

- .02876 

206.525 

206.506 

15 

.00000 

206.525 

206.506 

l6 

-.02313 

206.534 

206.507 

17 

.00000 

206.534 

206.507 

18 

-.02113 

206.542 

206.508 

19 

.00000 

206.542 

206.508 

20 

- .04124 

206.576 

206.511 

21 

.00000 

206.576 

206.511 

22 

- .03279 

206.599 

206.524 

23 

.00000 

206.599 

206.524 

24 

-.07421 

206.731 

206.548 

25 

.00000 

206.731 

206.548 


^Methods used in approximating S(|), see the appendix. 
The theoretical solution for this body is: 

00 

^ nAn^ = 2Aa^ = 206.521 sq. in. 
n=x 
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TABLE H.- TCHEBICHEF SOLUTION FOR AN AREA CURVE 
WITH A DISCONTINUITY OF SLOPE 


.4 


SCf) 




1.0 




n 

1 

II 

Tn 

-0 .6 
-1 -0 

An( machine computed)^ 

1 

- 0.254648 


- 0.254648 

2 

.407437 


.407438 

3 

-.410832 


-.410837 


.29335^ 


.293359 

5 

- .136980 


- .136983 

6 

.026293 


.026296 

7 

- .001958 


- .001961 

8 

.045998 


.045999 

9 

- .104130 


- .104131 

10 

.126592 


.126594 

11 

- .099202 


- .099202 

12 

.046071 


.046071 

13 

-.006268 


- .006269 

ll 

.003874 


.003874 

15 

-.032859 


- .0328^ 

16 

.065408 


.065407 

IT 

- .074838 


- .074839 

18 

.054970 


.054970 

19 

- .022337 


- .022337 

20 

.001456 


.001456 

21 

- .005707 


- .005707 

22 

.028363 


.028362 

23 

- .049488 


- .049490 

2h 

.052133 


.052133 

25 

-.034900 


- .034900 


1.0 


^Macliine computation, linear approximation of 
S(|)^ see the appendix. 




TABLE III.- COMPAEISOW OE TCEEBICHEF SOLOTIOES WITH VALUES COMPUTED BY 
THE METHOD OE EEFEEIENCE 1 EOR A BASIC SWEPT-WIHG MODEL (REF. 3) 


Area Curve 



m=i. 5. 

9 = 90 °, 1 = 0 ° 

M=l-5, 0=22.5°, ^=46.0° 


TcLebichef 

Method of ref. 1 

Tcheblchef 

Method of ref. 1 


N 

N 

N 

H 

N 



I 

Z 


n=i 

n=i 

n=x 

n=i 

1 

0 

0 

0 

0 

2 

282 

277 

278 

272 

3 

283 

277 

278 

272 

4 

315 

315 

284 

284 

5 

316 

317 

290 

289 

6 

405 

390 

291 

290 

7 

4l4 

396 

322 

320 

8 

521 

513 

326 

326 

9 

529 

524 

344 

349 

10 

628 

610 

385 

393 

11 

643 

624 

386 

393 

12 

686 

676 

483 

480 

13 

703 

692 

493 

496 

l4 

722 

708 

533 

552 

15 

729 

717 

588 

602 

l6 

729 

717 

598 

605 

17 

732 

721 

649 

663 

18 

736 

724 

660 

676 

19 

737 

725 

680 ■■ 

695 

20 

749 

736 

712 

726 

21 

750 

736 

713 

729 

22 

757 

744 

738 

75^ 

23 

759 

744 

759 

781 

24 

763 

748 

759 

781 

25 

764 


783 



Ho discontinuity of slope. 
Difference ~ 2 percent 


Discontinuity of slope. 
Difference 3 percent 


See figure 4 for plots of TcHebichef cbeck solutions. 
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TABLE IV.- FIBELAGE OEDINATES OF BASIC BODY - MODEL 3 


FiLselage station, 
in. 

Eadius, 

in. 

0 

0.28 

5.00 

1.46 

10.00 

2 . 4 i 

so. 00 

3.90 

30.00 

5.07 

40.00 

6.01 

50.00 

6.78 

60.00 

7.40 

70.00 

7.86 

80.00 

8.20 

90.00 

8.41 

102.00 

8.50 

110.00 

8.46 

120.00 

8.30 

130.00 

7.99 

140.00 

7.67 

150.00 

7.36 

160.00 

7.10 

170.00 

6.70 

180.00 

6.27 

190.00 

5.77 

200.00 

5.25 

210.00 

4.50 

219.00 

2.32 

225.38 

0 



36 




WACA EM A55J28 


TABLE V.- FUSELAGE ORDINATiS OF SCOOP BODY - MODEL 3 
[All dimensions given in Indies.] 


station 

Dlaenslons 


SHaezislons 

r 

1 y 

X 

I 

z 


r 

K 

X 

T 


25.00 

1353 

■■i 




111.31 

8 .t 5 


0 

2.00 

t.OO 

6.00 

7.00 

7.50 

7.75 

8.00 
8.25 

8.50 

8.95 

8.75 

15.90 

15.6a 

14.76 

13.12 

11.63 

10.67 

9-35 

7.58 

2,85 

0 

1 

ms ^ 

msm 

■Mi 





30.00 

5.07 


H 

B'J&] 

IB 


l^llllllll^ 

6.78 


1.50 

3.00 

t.76 

6.00 
6.75 

6.78 

5.85 

5-72 

t.97 

3 .t 3 

.82 

0 

1 

70.00 


1 

0 

3.00 

5.50 

6.50 

7.50 
7.80 

8.00 
8.05 
8J.7 
8.23 

^.32 

6.32 

6.32 

6.23 

5.62 

tTS 

0 

1.25 

2.50 

i.ti 

.39 

120,88 

8.28 


0 

2.00 

t.OO 

6.00 

7.00 

7.38 
7.66 

8.00 
8J.5 
8.25 
8 .to 

8.39 

15.90 
15.62 
lt.76 

12.91 
10.61 

8.10 

6.00 

3.98 

2.18 

.33 

0 

3.50 

2.18 

1.57 

1.15 

1 

&.ao 

it.53 
It .17 
13.18 
U.21 
9 . 5 ^^ 

8.39 

7.70 

6.to 

0 

2.00 

t.OO 

6,00 

7.00 

7.50 
7.75 

8.00 
8.25 

8.50 
sM 
8.7? 
8 . to 

6.^ 

6.52 

6.52 

6't8 

6,39 

6.09 

3 .t 6 

2.03 

0 

0.80 

130.00 

6.02 


0 

2.00 

t.OO 

5.80 

6.25 

6,30 

6.50 

7.00 

7.50 

7.75 

8.00 
8.05 

15.90 
15.62 
14,76 
12.98 
10.50 
7.30 
5. to 
t.03 

ee<*e» 

2.20 

.95 

0 

2.82 

2.07 

.79 

83.00 

1 


0 

1.50 

2.50 

3.50 
t.50 

5.00 

6.00 

7.00 

8.00 
8.25 

8.50 

9.00 
8.75 

8.57 

15,15 

14.98 
lt.68 
It .at 
13.60 
13.19 

12,17 

10.80 

8.80 

.92 

0 

2.25 

l.t 5 

.45 



i 5.90 

15.62 
It, 76 
13.76 
13.18 
12.81 
12.00 
6.70 
6.00 
5,26 

0 

2.00 

t.OO 

5.25 

5.60 

5.70 

5.80 

5.80 

5.85 

6.00 

5.00 
5.50 

5.85 

6.59 

6.11 

5.96 

5.98 

5.96 

5.90 

5.90 


102.00 

1 

1 

0 

2.00 

t.OO 

6.00 

7.00 

8.00 
6.25 

8*^ 

9.00 

9.35 

9.00 

8.97 

15.90 
15.62 
lt.76 
13.12 
11 .8t 
10.10 

8. go 

7.20 

3.70 

•15 

0 

2 .tt 

ito.oo 

7. 60 


0 

1.92 

3.74 

4.28 

t .88 

5.37 

1 

1 

jj^m 

1 

1 

1 

1 

1 


^ i ^ » ■ mi I ■ I . , , 4 

Dote: See figure 11(e) for definition of sj^mbola. 

See tstle 17 for redU of nose and tall portlona of fuselage. 
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Figure 1.- Definition of the Tchebichef polynomlalB used to represent the Fourier sine series 
defining the slope of the area diatrlbutlon curves of a model. 
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IS •['-(f)']' 

S'(s(} = Aesin2(^ 


CdJ= f 

do “ 2fo 


Area 

distribution 



X — * 

'[ w 

All dimensions are in Inches 

ro = 

26.96 r 

^ . X 


Body shape 

Figure 2.- Ml 3 oinmni-drag body for preacrlbed voliinie and length (Sears-Biaack body); ftneneea 

ratio 12. 5. 
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Figure 4 




IW 


Fuselage station, in. 


(a) Kb discontimity of slope, M = 1*5; 9 = t = 0^- 
- Plot of Tclxebiclief elieck solutions for area curres previously analyzed (ref. 3 ). 





Fuselage staNon, in. 


(b) DlBcoirfcliniity of slope, M = 1.5^ 0 = 22,5°^ if 


mi<na 1 l ^ 
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Figure 5 .- Trlangular-wiag intercepbor-type airplane model (ref. 9)j Model 1. 
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Experiment, ref. 9 

Computed, wave-drag coefficients, N<»25 
Computed, friction-drag coefficients 
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0 !0 20 30 40 50 60 70 80 90 !00 120 

Fuselage stotion. In. 

Figure 6 .- ^Eriangular-wing airplane (JfocLel 1, fig. 5 ) tnlnitnum-drag coefficiente aaad cross- 

sectional area distribution. 





Cross-sectional area, S^, sq in. 



(a) M = 1.00 

Figure 7.- Projected-area distribution for the triangular-wing airplane (Model 1) for Mach nujnberB 

of 1.0, 1.2, and 1.8 with check nolutiona for N = 25 . 
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Fuselage station, in. 


o 


(b) M = 1.2; 0 = 90°; t = 0° 
Figure 7.- Continued. 
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Fuselage station, in. 

(d) M = 1.2, e = 270°, + = 0° 

Figure 7*" Continued. 
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aireOi S, sq in 




SVa 9.6 




WtNG 

Aspect ratio, 3.56 
Area fiofaii, 376,02 sq ft 
Dihedral angle. 0 deg 
M.AD. 11.33 ft 

Airfoil section, NACA 64A007 

HORIZONTAL TAIL 
Aspect ratio, 3.56 

OQ ere 

Miow |iw«wir| «««| ■■ 

Dihedral angle, 0 deg 
iyU.C, 5.81 ft 

Airfoil section, NACA 64A007 

VERTICALTAIL 
Aspect ratio, 1.27 
Area (BKCluding dorsal fin), 
48.35 sq ft 
M.AC. 6.61 ft 

Airfoil section, NACA 64A007 



All dimensions in inches. 

All airfoil sections parallel to the 

linA 

vvsiiwi <i«iw* 
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Figure 8.- 


Swept-wing Intexceptor-type airplane similar to low-tall version of reference llj 

Un>l.n-| O 
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Fuselage station, in. 



Ml 

U> 


Figure 9*"’ Croas-sectional area distributions of low-tail versions of a swept-wing airplane 

(O.U scale models). 
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Zero-lift drag-rise coefficient, Zi Co 


Experiment, low-tail version, ref. II 

0 Computed, low-toil version, ref. M 
El □ Computed^ airplane model 2 





Mach numoer, M 


Figure 10.- Comparison of exper imant Al drag- rise coefflclen'ts frcin reference H with compirted 
wave-drag coefficients (]J = 25 ) for two low- tailed versions of a swept-wlng airplane. 





Body radii are given in table IV 
External wetted orea IIOI6 in.^ 



Section AA (typical) 


(a) Basic body and tail. 

Figure 11.- Scoop-inlet duct model; Model 3 . 
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Figure UL.- Continued. 
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External zero-lift drag coefficient, Zero-lift drag coefficient 
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Mach number, M 

(a) Total and internal drag coefficient. 



.80 .84 .88 .92 .96 1.00 1.04 1.08 1.12 1.16 

Mach number, M 

(b) External drag-rise coefficient. 


Figure 12.- Experimental and computed drag coefficients for the scoop- 

Inlet-duct model; Model 3 . 
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Fusetage station, in. 


(a) Inadmissilole check solutiona for K = 10 and N i= 25 . 

Figure 14.- An Inadmissible area curve with infinite slopes. Scoop-inlet duct model for duct 

mass-flow ratio of 0.6j Model 3 . 












